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When using black hole excision to numerically evolve a fully generic black hole 
spacctime, most 3-D 3 + 1 codes use an :rj/2-topology (spatial) grid. In such a grid, 
an r = constant excision surface must be approximated by an irregular and non- 
smooth "staircase-shaped" excision grid boundary, which may introduce numerical 
instabilities into the evolution. In this paper I describe an alternate scheme, which 
uses multiple grid patches, each with topology {r X (angular coordinates)}, to cover 
the slice outside the r = constant excision surface. The excision grid boundary is 
now smooth, so the evolution should be less prone to instabilities. With 4th order 
finite differencing, this code evolves Kerr initial data to ~ 60M using the ADM 
equations; I'm currently implementing the BSSN equations in it in the hope that 
this will improve the stability. 

1 Introduction 

Black hole excision (BHE)§! is now widely used for the 3 + 1 numerical evolution of 
dynamic black hole spacetimes. In 3-D, most codes use rryz-topology (spatial) grids 
to avoid z axis problems, with a "staircase-shaped" excision surface approximating 
r = constant. Unfortunately, the irregular shape of such an inner boundary to the 
grid considerably .complicates the finite differencing, and often introduces numerical 
instabilities (e.g. 0). 

2 Multiple Grid Patches 

To avoid the numerical problems caused by irregular grid boundaries, I pro- 
pose a new scheme, in which multiple grid patches, each with topology {r x 
(angular coordinates)}, are used to cover the slice outside the r = constant ex- 
cision surface. This allows the inner boundary to be smooth. By using multiple 
patches, the angular coordinates can be non-singular in a neighborhood of each 
patch, yet still cover the full angular domain. 

Figure |l| shows an example of a 3-patch coordinate/grid system of this type, 
covering one octant of space. There are many possible ways of constructing such a 
multi-patch coordinate/grid system; the particular choice used here has the advan- 
tage that each pair of adjacent patches shares both r and and one of the angular 
coordinates. Each patch is surrounded by angular "ghost zones" (not shown in the 
figure); field variables in the ghost zones are interpolated from neighboring patches. 
Because of the shared coordinates between adjacent patches, only 1-dimensional 
(tensor) interpolation is needed. By overlapping the patches a few grid spacings, 
the interpolations always use data from patch interiors, which helps the evolution's 
numerical stability. 
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3 Numerical Results 



I have implemented this scheme in a new BHE evolution code (about 35K lines 
of Maple and C++), using 4th order finite differencing in space, a smoothly- 
nonuniform radial grid, and a 4th order Runge-Kutta method of lines time in- 
tegration, all basically as in my spherically symmetric code (□). With the standard 
ADM equations, the code currently has finite differencing instabilities at the corners 
where all 3 patches meet, with the energy constraint reaching ~ f % by t ~ 60M fox 
Kerr initial data with spin 0.6. I'm currently implementing the BSSN formulationta 
and investigating changes to the interpolation scheme to try to improve the stability. 
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Figure 1. This figure shows an r = constant shell of a 3-patch coordinate/grid system covering 
one octant of space outside the excision surface r = r m - ln . The patches overlap by 3 grid spacings. 
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